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Vor. III. OcToseR, 1887. No. 5. 


ON THE CONDITIONED COOLING AND THE CUBICAL CONTRACTION 
OF A HOMOGENEOUS SPHERE.* 


By Mr. R. S. Woopwarp, Washington, D. C. 


13. In the former paper (ANNALS oF Matuemarics, Vol. IIT. No. 3) we have 
discussed the problem of the “free cooling” of a homogeneous sphere. In the 
present paper it is proposed to consider what we have, for brevity, called the case 
of “conditioned cooling of a homogeneous sphere,” or the cooling of such a 
sphere in a diathermous medium. In addition, we shall give an application of 
the theory of cooling in determining under certain limitations the consequent 
cubical contraction in large spheres, this problem being of special interest with 
reference to the secular contraction of the earth. 

The physical circumstances, and the mathematical features of the solution 
of the problem of conditioned cooling have been very fully and ably set forth by 
Fourier,t} Poisson,{ and Riemann.$ So far as the purely analytical requirements 
are concerned, their works, especially that of Poisson, will be found exhaustive. 
We need, therefore, only to refer to them for a complete exposition of the ele- 
mentary considerations which lead to their final result, our object being to adapt 
that result to the requirements of the computer, having in mind always an ap- 
plication of the theory to the earth. It will be desirable, however, for the sake 
of clearness and definiteness to make a formal statement of the fundamental as- 
sumptions and relations involved in the steps of the investigation which precede 
those we have in view. 


* Since this paper is essentially a continuation of the one on “ The Free Cooling of a Sphere,” pub- 
lished in Vol. III. No. 3 of the ANNALS OF MATHEMATICS, the section numbers and those designating the 
principal equations are made continuous with the corresponding numbers of the earlier paper. The nota- 
tion introduced in the latter has been adhered to as far as possible. 

+ Théorie Analytique de la Chaleur. Paris, 1822. Analytical Theory of Heat, by Joseph Fourier. 
Translated, with notes, by Alexander Freeman. Cambridge, 1878. 


t Théorie Mathématique de la Chaleur. Paris, 1835. 


2 Partielle Differentialgleichungen und deren Anwendung auf physikalische Fragen. Vorlesungen 
von Bernhard Riemann, fiir den Druck bearbeitet und herausgegeben von Karl Hattendorff. Dritte Auf- 
lage. Braunschweig, 1882. 


ANNALS OF MATHEMATICS 
e 

j 
43 


130 WOODWARD. ON THE CONDITIONED COOLING 


14. The assumptions made are four, viz.: first, a constant coefficient of dif- 
fusion for the entire mass of the sphere; second, an initial uniform excess in tem- 
perature above that of the surrounding medium of the entire mass of the sphere ; 
third, a rate of flow of heat from the surface of the sphere to the medium pro- 
portional to their difference in temperature ; and fourth, a constant temperature 
for the surrounding medium. 

As in the previous paper, let 7 be the radius of any spherical shell whose 
thickness is dv, « the excess in temperature over that of the surrounding medium 
of this shell at the time ¢ after the initial epoch, and a’ the coefficient of diffusion. 
Then these quantities must satisfy the partial differential equation 

(ru) 2°? (ru) 
Ort 


(29) 


cr 


and the solution of the problem consists in determining the integral of this equa- 
tion, (7) being regarded as the unknown quantity, subject to certain relations 
which we proceed to set down. If be the initial uniform excess of the tem- 
perature of the sphere over that of the surrounding medium, we must have 


ru = ru, for = 0. (30) 
In order that ~ may not be infinite when y = 0, we must have 
ru = 0 for r =o. (31) 


It remains to express the characteristic condition of the problem, namely, 
that which defines the emission of heat at the surface of the sphere. If & be 
the coefficient of conductivity for the sphere, the quantity of heat which passes 
across a unit of its surface in an element of time is expressed by the value of 


—k x when + = 7,, the radius of the surface of the sphere. On the other hand, 


this quantity of heat is equal to that which the surrounding medium can carry 
away, or according to our third and fourth assumptions, to Hu, * being a con- 
stant and the temperature corresponding to r= Hence we have 


ou 


—ky = Mu, 


or, if we put 7/4 


A 
cu 
+ hu=o forr=r). 


(vw) Ou 


But since = 


*This quantity is cailed by Sir W. Thomson the emissivity. It is the quantity of heat emitted per 
unit of time, per unit of surface, per unit of excess of temperature; and may be due to conduction, con- 
vection, or to both. See 771, Article Heat, Encyclopedia Britannica, 


i | 
a 
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which, substituted in the above relation, gives 


(ru) 
ov 


+ (rh —1)u =O forr=7,. (32) 


Now a particular integral of (29) satisfying the relation (31) is 
ru = C'e-*"™ sin Ar, (33) 
in which C’ is a constant, and / a factor such that 


(ru) 


= cosir 


and ua Clr-'e-*"™ sin dr 


will satisfy (32). Substituting these values in (32) we have 4 defined by the fol- 
; lowing transcendental equation : — 


r,Acos rh + (hk — 1) sinr,A =. (34) 


This equation has an infinite number of unequal real roots,* and since any one 
of them will satisfy (33), the complete integral of (29) must consist of a series of 
terms of the form (33). That is, 


ru sin dy. (35) 


In this, A,, 4, 4s, ... are the successive roots of (34). The constants C’,, C’,, 
C’,, ... must be such that (35) will satisfy the initial relation (30); i. e. making 
¢= 0 in (35), the relation 


© 


The process of determining C’, from this equation consists in multiplying both 
sides by sin A, 7dr and integrating between the limits 0 and 7, account being 
taken of the relation (34). We shall not go through the steps of this process, 
referring the reader to Poisson and Riemann for the details, but shall pass imme- 
diately to one of the usual final forms for (35). This form ist 


0 
2u 

ru = 73 
A,, 


2) 2 2 

) ga atta n? (sin r,A,, cos sin Th, 
+ (rh — 1) 


(36) 
15. The expression (36) satisfies all of the relations (29) to (32), because each 
of its terms satisfies them, and is the complete solution of the problem. But, 
though it meets all the demands of pure analysis, its form is quite repulsive to 


* See any of the works first cited above. 
+See Riemann, Differentialgleichungen, 769. ther equivalent forms are given by Fourier, Poisson, 
and Riemann, 
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the computer or to any one who may desire to use it, considering that the quan- 
tities 4, must be determined from the transcendental equation (34). Some at- 
tempts have been made to simplify this expression, but the results are not satisfac- 
tory, inasmuch as they require either a time so great that we may neglect all terms 
of the series after the first, or a very small radius of the sphere, or a very great 
radius of the sphere. No attempt, so far as we know, has been made to trans- 
form the series into another more converging or tractable, though Poisson has 
summed it to a limited degree of approximation, as we shall presently explain. 
Neither is the improvement we have to suggest entirely satisfactory, for it is re- 
stricted to the case where 7,/ of (34) is greater than unity,* but this is the case 
with which we are most particularly concerned and that to which we now direct 
attention. 

16. The quantity /% in (34) ts of — 1 dimension in length,t and hence 7,/ is 
of zero dimension or an abstract number. 

For brevity, let us put 


Then equation (34) gives 
tan y, = (37) 


Likewise, if we make the same substitutions in (36), it may be put in the more 


compact form 
= 27, 72° + xr 


Now it will be observed that when + = 0, v7, = uz, 4, = nz/ 7, and 


sin (38) 


siny, 
a 
hence (38) becomes 
27,4, « (—1)*t! r 


This agrees with equation (10), which expresses the solution in the case of 
free cooling. The case of conditioned cooling, therefore, approximates to that 
of free cooling in proportion to the smallness of the number x This number 
may be small by reason of the largeness of either 7, or 2, or of both. From cer- 

*It is obvious, however, that a slight modification of the process followed in the text may be used in 


treating the case where 7, 4 <1. 
¢ See Fourier’s Section IX. chap. II. p. 130, Freeman’s translation. Our / is the same as Fourier’s 


h| K of the page just cited, but the same as his / defined on p. 268. 
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tain observations of underground temperatures in the garden of the Paris Obser- 
vatory, Poisson concluded that for the earth in the vicinity named, / is a little 
greater than I, using the metre as space unit.* This would make x = 1 (”,/) 
about 1 / 6,300,000. Whether this value is nearly correct as applied to the whole 
surface of the earth may be doubted, but + might be 10,000 times as great with- 
out materially affecting the applicability to the earth of the process we shall 
develop. The correctness of that process only requires that + < 1. 

It is proposed, now, to expand the second member of (38) by Maclaurin’s 
theorem in a series of ascending powers of +. As a preliminary, though not es- 
sential part of this work, we may expand 1, of (37) as a function of x and arrive 
thus at a formula of use in a future numerical application. 

17. When + = 0, (37) gives 7, = wz, and this is the first term in the expan- 


a‘ 
sionof y,. To get the coefficients of the terms involving powers of 1, we need ra 


for + =0,¢ being any positive integer. If we differentiate (37) once, and then 


eliminate sec? y,, the result is 
(1 + +), =0 


If now we put Y= y, : ", and differentiate the above equation successively, we 


find that for + =o 
4 
6-9 Sr (39) 


By a duplicate application of the theorem of Leibnitz, we can derive a gene- 
ral expression for the (2 — 3)th differential coefficient of y with respect to +, but 
its usefulness hardly justifies the space required by its rather complex character. 
We shall simply write down the first six values of (39). They are 


ay, 


nz fort= 


= 
=-+ 
=— 8 
=— 
= + 384 (uz) 
Hence by Maclaurin’s theorem we have 
In = (1 — + — — 


*Théorie de la Chaleur, p. 502. The value given is 6= 4 = 1.05719. 
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Neglecting terms involving powers of « higher than the 4th, this equation may 
be written 

Vy = (A — (1 + 

= nz — {1 — (1 — 2)}, (41) 

one of which forms will be found convenient for logarithmic computation of the 
smaller values of y,. For the larger values, the series (40) will, in general, prove 
too slowly converging, but it is not necessary to calculate such values. Our only 
use for the smaller values is a numerical verification. 

18. To develop (38) we observe that each of the seven factors under the sign 
of summation, is a function of x. We might, therefore, expand these factors sep- 
arately and then find their product. On the other hand, we might get the differ- 
ential coefficients of (77) with respect to 2 by direct differentiation. We shall 
adopt a combination of these two processes. 

First, develop the product +~'sin_y,. By means of (40), we readily find 


a~'sin y, = — uz cos uz (1 — + 
Then, if we put y,, = xzy, where y is the series within the parentheses of (40), we 
see that 


—1)*t! |. 
sin y, yr — + 


Now, for brevity, let 
+ 
and X, = (1 — + 
Then, if we denote the ~” term of the series under the sign of summation in (38) 
by 7,, and differentiate log 7,, the result is 


T, ox Y ade a. 4; ax +s ax 
ah 
2a‘th,, + rcots,r 
dy 
But from (40), we have when x = 0, y = 1, and 3, 8. Also, since 4, = nf 
0 
ux dh, 
and 
Therefore, when x = 0, we find 
= 7, | 1 + 2¢ cot uz |. (42) 


The succeeding complete values of the differential coefficients grow rapidly in 
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complexity. Their derivation, however, is tedious rather than difficult. It will 


suffice, for our purposes, to write down the special value of the second differen - 


r? anz)* ) 
r 


tial coefficient, viz. : 


r nz) r 
2—ar+ 4alr| — cot — 


Now we observe that when + = o, 


( 1)" attint (ro)? os 
NX 0 
* 


r é 
7, — nx 
a 


(¢) 


For brevity, put S for the value of the second member of (38) when + = 0, this 
value being the same as that given by (10), (19), or (20). Then, observing the 
effect of the sign of summation in (38), equations (42) and (43) in connection 


with the equivalents in (a) to (/) give 
(ru os Os 
(44) 
oS 
4S — 2r 


+ 


*To get the signs of the second members of (4) and (/), one must note that 97, / 97 and da, | dx 


have opposite signs. 


) 
(43 
& 
(a) 
| 
4t? | 
227° | | 
2 
(d) 
’ | 
2 
roe > 
0 
+ 
I 4 4 a 
re: ( 
& 
Ta 
4 
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Finally, by Maclaurin’s series we get 
ru = S+ (46) 


19. The partial differential coefficients involved in the second members of 
(44) and (45) are easily calculated. If we write, for brevity, 


( gh — ity — we)? ) 

(47) 
n=I1 + — + +m) 


n=0 ( (2m — m,— m) 

= + m, + m) Cmontmotm) 
( (2m — m, — 

Oo!” 

= 2 + +m)? (49) 
( + + my, + e~ Amon timo 
(2m — m,— my ) 

Oov= -(50) 

— (2iton + + m)2 (5 ) 

— (21,2 + m, + mye 


oS 
=+r, (1-5 0'), (51) 
Os 
2t = — (52) 
2 Vo? i 
= + — 40"), (54) 
— + (20’ — 40”). 
+ 4 (55) 


These will be applicable for all but the latest stages of cooling, and the corres- 
ponding derivatives from (10) will become applicable when these fail to be suffi- 
ciently converging. 

The solution for conditioned cooling is thus made to depend on that for the 
ideal case of free cooling; and the expression (38), which is practically useless 


m, and m being defined by equations (18), we get from (23), (24), and (25) 
| 
J 
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(on account of its slow convergence) for smal! values of the time is transformed 
into (46), which is available for any value of the time. We have limited the ex- 
pansion to terms of the second order in x, believing that this is more than suffi- 
cient for any application of the theory to the earth; but the higher terms of the 
series may be easily derived if needed. 

20. To complete the solution, it remains to show how to determine the tem- 
perature at the centre of the sphere. When we make ¢ = 0, the expression for 


w from (46) assumes the ambiguous form * but the usual process of differentia- 


tion gives 


af (v7) (rv) 
oS or ox* 
(56) 
er er cr 
From equation (26) we have 
a” 
S m 
cr 
r= 
The differential of (44) gives 
OF eS e*S 
= — (58) 
cr or crete 
r =O. 
From (55) we have 
rou 
eret ay (59 


in which, of course, ” must be made zero for the special value in (58). 

The special value of the differential of (45) with respect to , required by 
the second term of (56), may be easily written down and worked out numerically 
for any assumed data, but the paucity of practical applications makes this work 
unnecessary. * 

*In note C of the supplement to his great work (see full title, p. 129) Poisson has summed the series 
(38) to terms of the first order in rv. He does not separate his result into the term independent of and 
that dependent on .r, but gives an expression for the sum of the two terms. He deduces several special 
results dependent on the assumption (@) that x = 0, or (4) that the strata whose temperatures are consid- 
ered are near the surface of the sphere. Among these results is an expression for the temperature at the 
centre of the sphere under the assumption that «0. This expression agrees with our expression (26) or 
(57); but the casual reader would not observe the agreement; and, curiously enough, Poisson himself was 
unaware that the difference of two definite integrals in his expression is exactly zero. Using our notation, 
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21. Since no set of formulae presenting any degree of complexity, or requir- 
ing many steps in their derivation, can be regarded as entirely trustworthy until 
they have been subjected to the test of actual application, we propose in the pres- 
ent and next sections to work out in detail a numerical example. For this pur- 
pose we shall assume such a value of the time that we can compute +z from (38) 
as well as from (46). The equivalence of the two forms will thus be tested, as 
well as the correctness of the signs of the several terms in the second members 
of (44) and (45). 

We shall assume m, = 7, / 2a,/¢ = 7/2, a value used also in the illustrative 
example of $11. This gives 


which, for the earth, if we use Sir W. Thomson's value of a’, is, as stated in $3, 
about 100,000,000,000 years, With this value of the time, the first three terms 
of (38) will suffice to give vu in terms of 7,7, to seven places of decimals. We 
shall also assume + = 10~* and y = $y, and compute first by means of (38). 

22. To begin with, we need the first three roots of the transcendental equa- 
tion (37). These we get from the first of (41), which gives with sufficient accu- 
racy log y,, = log (wz) + log (1 — +) + 4 (uz) x*; » being the modulus of the 
common system of logarithms. This gives 


log y, = 0.4967154, 
log 7’, = 9.797745, 
log 7’; = 0.973837. 
For use in sin y, and sin 4,7, it is best to compute the supplements (7z — In): 


he has for the fall of temperature at the centre of the sphere (see p. 52 of the supplement) 


“uy 0 [/ dz + 2m, — Mo? — 1)? — + =)? dz| : 
Mo n=l n=1 —mo 
while our equation (26) gives only the middle term of this expression, But if in the last definite integral 
we replace 2,2 + z by 2, the result is 


ity (2n+1) 3mo 


1 my (2n — 1) mo 


= ds. 
me 


which is identical with the first definite integral. Hence the above expression is the same as (26). Pois- 
son, after reaching his equation, observes that 


rn 


+ 
— mon + dz — (2H — 


and proceeds to draw some conclusions which are qualitatively correct, but not susceptible of definite nu- 
merical evaluation. 


t = _v 
ar 
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These may be accurately computed from the second of (41), which gives in sec- 
onds of arc 
— = 647797787, 
27 — J, = 1295.983, 
37% — = 1943-94- 
These and the above values of log y, will be found to satisfy (37). 


We can now compute all the factors in (38). The logarithms of the second 
algebraic factor and those of the factor x~' sin y, are found to be 


log log sin y, 
9.9995055 — 10, 0.4971479, for x = 1, 
9.999505 — 10, 0.798170, == 2, 
9.999506 — 10, 0.974252, == 3. 


Since 4, = y,,/7%), we have with our assumed value of / 


Ym. 7)? 


and hence log = 9.56657 36 — 10, for 7 = 1, 
= 8.266294 — 10, = 2, 
= 6.099162 — 10, = 3. 
Likewise, for the last factor we have sin 4,7 = sin y, 7/7, = sin }y,. Hence 
log sin }y, = 9.9999995 — 10, for x = 1, 
= 7.798170 — 10, = 
= 9.999995, — 10, 
Summing the logarithms of the factors of the several terms we get 


bo 


log No. 
Ist term, 9.0698557 — 10, + 0.11745073, 
2d 4.905681 —10, + . 0924, 
3d 5.125300, —10, — . 1334. 


The sum of the terms is 0.11744663. Multiplying this by 4 we get as the result 
from (38), 
= 0.4697865%,. 
23. The details of the computation by equation (46) are as follows : — 
The value of twice the first term, S, has already been computed in the ex- 
ample of §11; i. e. 
S = 0.23417 327%) 
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using the larger of the two values there given, since a more extended computa- 
tion shows that to be the more precise value, 

We next need the values of the Q's which enter the second members of (51) 
to (55). Since by (18), # =) 2a, ¢, and since for this example we have assumed 


Mm, 2a, and = we have for use in (47) to (50) 
Mm m+m= Zz, 
3m, m= 3r, fz, 
5m, — m= 5m, + m= iz, 


Hence we have, to terms of an order which may be neglected, the following sys- 
tem of equations in which the exponential quantities are factors respectively of 
the several coefficients in the vertical rows: 


== + Gay + (fey + 


These give, by an easy computation, 
QO = + 054352, 
Q” = + 0.41469, 
ON” = + 0.35443, 
OV = + 0.21069. 
Introducing these factors and the values of ¢ and y in equations (51) to (55), we 


find the following equivalents, in which 7,7, is a common factor of the second 


members : — * 


r = + 0.018327,%,, 
~ 
2t a, = — 0.46793, 
(r2 = — 2.886, 
PEAY 
40? = 0.928, 
PEAY 
4rt Dror = — 0.293. 


*The numerical values of all these differential coefficients have been checked by means of the equiv- 
alent forms obtained from equation (10). 
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By means of the last values and the above value of S we get from (44) 
© (ri) 
Likewise from (45) 


-3 


With these values and with + = 10 *, we have from (46) 


u == 2u, (0.2341732 + 0,00072042 — 0.00000044) 


= 0.46978647,,. 
This agrees with the value obtained from (38) as closely as 7-place tables of log- 
arithms will permit, the discrepancy being one unit in the last place. 

It is scarcely necessary to remark in this connection that the precision of the 
preceding computations 1s of little importance except as a test of the correct- 
ness of the formule derived and as an illustration of the process which might 
be followed if our data for practical applications were sufficiently trustworthy. 
For all applications to the earth, computation of results by the preceding for- 
mula may be properly limited, in general, to three or four places of decimals, 
since defects in the best data at present attainable and probable inaccuracies in 
our fundamental assumptions (of $14), render those results uncertain even in the 
first place of decimals. For these reasons, also, and in view of the probability 
that the quantity + is less than 1/ 1000, it is believed that the solution for free 
cooling is an ample approximation in all applications to the earth. 

24. In this and the following sections it is proposed to determine the amount 
of cubical contraction of a homogeneous sphere on the assumption of free cool- 
ing, and on the additional assumption that the coefficient of contraction is con- 
stant. These assumptions will give correctly the first approximation (or the term 
independent of +) in the case of a homogeneous sphere of moderately high initial 
temperature, and probably as good an approximation as is at present attainable 
for the secular contraction of the earth.* 

If ¢ denote the coefficient of cubical contraction and J!’ the total change in 
volume during the time ¢ after the initial epoch, we shall have 


ro 


4V= are (% — u) rdr, (60) 


in which the symbols under the sign of integration have the same meaning as 


*The consideration of the mechanical effects of secular contraction forms no part of our present 
work, though one of the principal objects of the analysis of this and the preceding paper is to furnish a 
basis for the investigation of those effects. 


4 
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heretofore. The fall in temperature, 7”, — #, is given by equation (21) or (22), 
but for our present purposes it will be convenient to use another expression easily 
derived from (10). If in the latter, we make ¢ first zero, and then ¢, and take the 
difference between the two expressions the result is 


0 
| 


Introducing this value of the fall in temperature of the shell considered, in (60) 


we get 


JIV=8r,u,¢ (1 e rdr sin nz 
n 
Oo 
But the integral in this equation is equal to 
_%y COSNT ( 


Therefore we have as the exact solution, under the assumed conditions, 


When ¢= 0, this expression gives JV = 0, as it should. When ¢= «, the 
quantity under the sign of summation becomes 


34+... .=7'/6; 
and hence IV = 


which is also obviously correct. 

25. For large values of the time the series in (62) is evidently rapidly con- 
verging, but for small values the convergence is too slow for practical applica- 
tions. It is desirable, therefore, to transform (62) and derive an expression from 
which we can compute the contraction when ¢ is small, or during the earlier stages 
of cooling. To this end we proceed. 

Referring again to the Gamma function of §4, and making the same substi- 
tutions for « and # as in §5, the series in (62) becomes 


4x? — a dx (cos zx + 4 cos 27x + +. ..), 


in which x is merely a subject of integration and is, of course, not the x of 


AND THE CUBICAL CONTRACTION OF A HOMOGENEOUS SPHERE. 143 


$$16 to 23. Now, for brevity, denote the trigonometric series in the last expres- 
sion by a,, and let 

do, a*a, da, 

ax 


From the series in $6 we see that 


o = — ( When g = 1. 
I — + ¢ 


Multiplying by 2, integrating between the limits o and x, and making g = 1, we 


get 
6, = 47 2. 


Multiplying the last equation by @v and integrating again between o and 4, there 
results 


a, = (§ + — fa). 
This gives the value of @, for all values of x from o to 2 both inclusive, or through 
the first period. For succeeding periods we have 
[4 + 2)? —4(4 — 2)], for x between 2 and 4, 
= + 4 (4 — 4) —4 (4 —4)], for x between 4 and 6, 


Substituting these values of 4, in the above integral, and writing, for brevity, 
m, = %,/2a,/t, as in equation (18), we get 

} 


— 2m,73 


4 
+ +4 (4 — —4 (4 — 2)] 


6 
4 


Collecting and arranging the terms of this expression and performing the inte- 
grations so far as possible, it reduces to 


an + 2 
0 
2n 


Substituting this value for the quantity under the sign of summation in (62) and 
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replacing x by m,, we get finally 


2m, 8, 
SV = - (63) 
+ [ fe “ds | 
My 2a) 


26. The series in the right-hand member of (63) converge with extreme sud- 
denness for small values of the time. For the case of the earth they converge 
so rapidly that their first terms alone will suffice for the most accurate computa- 
tions (one part in a million, say) when the time is as great as 7,?/ 4a’, which, with 
Sir W. Thomson's value of a’, is about 273,000,000,000 years. Using this value 
of the time we have computed the amount of the contraction by both (62) and 
(63). The first three terms of (62) will suffice in this case for the precision at- 
tainable with 7-place logarithms. The second term in the first series of (63) or 


e ' 4 x amounts to but six units in the eighth place of decimals; and the sec- 
ond term of the second series, or - 
6 
12 
fe dz, 


amounts to but twenty-eight units in the eighth place. Omitting other details of 
the computation, the results are 
JV = 0.1580729. from (62), 


= 0.1580728 . (63). 


These results serve to verify the numerical equivalence of (62) and (63), and show 
that they, like (46), or (10) and the group (19) to (22), bridge completely the 
time interval of cooling and enable the computer to give a definite answer to any 
question falling within the scope of the assumed conditions. 

It is worthy of remark that the contraction just computed is about 

SV = 0.95 . 427,50, ; 

that is, in the time ¢ =~," / 4a’, a sphere will, under the assumed circumstances, 
attain ninety five per cent. of its possible contraction. The temperature at the 
centre of the sphere at the end of this time is (see $12) about one-sixth its initial 
value. 

In all practical applications to the earth, for any but the latest stages of cool- 
ing, we may neglect the series in (63), which becomes then simply 


3 


From this it appears that in the earlier stages of cooling the contraction is pro- 
portional to the square root of the time. 


i) 


HILL. ON DIFFERENTIAL EQUATIONS WITH PERIODIC INTEGRALS, 145 


ON DIFFERENTIAL EQUATIONS WITH PERIODIC INTEGRALS. 
By Dr. G. W. Hitt, Washington, D. C. 


The independent variable being conceived as time, a system of differential 
equations may be said to admit periodic integrals when the values of the depen- 
dent variables either exactly, or with approximate tendency, after a certain lapse 
of time, repeat their series of values. In the latter case the larger the lapse is 
made the more nearly is the repetition brought about. Strange as it may seem, 
this subject, except in the case of simply periodic integrals, is, at present, not 
completely understood. The text-books on differential equations are almost 
wholly engaged with the cases in which, by certain artifices, the integration can 
be accomplished in finite terms or reduced to quadratures. In the treatment of 
physical problems, however, equations of this sort are rarely met with. Far 
more frequently it is found that methods of approximation must be resorted to. 
Cauchy appears to be the author who has done most for the elucidation of this 
part of the subject. His memoirs are in his later Zverceses and in the volumes 
of the Comptes Rendus for 1856 and 1857. In this article | propose to show how 
simply periodic integrals arise and afterwards to illustrate the general theory by 
treating a problem relating to the motion of a system of points. 


Having the independent variable /, and the two dependent variables + and .4, 
let us suppose the latter satisfy the equations 


dt 


A cross multiplication between the members of these equations gives 


The integral of this is, C being the arbitrary constant, 
= + C. 


The values of + and *, being known for a given value of ¢, we readily find the 
value of C proper to the special case we treat. By substituting the value of 2, 
derived from this equation in the first of the differential equations we get 


ars c]. 
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The expression under the radical sign is a function of 1; calling it , let us con- 
sider the equation .Y =o. Since real values of x are supposed to correspond 
to all values of 7, .\ can never be negative; and, from the way the constant C was 
determined, it is plain that, for the given value of ¢, Y is positive. Then in VY =o, 
let + be supposed to increase until a value 1 = dé is reached for which .V = 0, 
that is to say a real root of this equation. Similarly let + diminish from the same 
point until a value # = c is reached for which again XV = 0, that is a second real 
root Then, X being positive for all values of x which lie between ¢ and 4, if the 
latter are non-multiple roots, .\ is negative for values of + which lie just outside 
these limits. Thus « must necessarily remain within the limits c and 4. Also, 
in its motion, it always attains them; for suppose + is augmenting, then the rad- 
ical, which forms the value of dv dt, must be taken positively, and, from the law 
of continuity, must continue to be so taken until it becomes zero, that is until x ar- 
rives at the value 6. But @v/ d¢ cannot be positive beyond this point, for x can- 
not surpass 4. Hence, after this, the radical must receive the negative sign, and, 
consequently, « begins to diminish. Again, from the law of continuity, this 


‘diminution is kept up until + has arrived at the valuec. At this point the dimi- 


nution must change into an augmentation, for + cannot fall below c. Thus the 
movement of + is a continuous swinging back and forth between the limits ¢ and 4. 


We can put X= 


X being a function of «+ which remains constantly positive and finite for all values 
of + between ¢ and 6. We can then write 


RR 
dx ~ 


A new variable ~ can now be advantageously introduced in place of x Let 

4+ =a(I —ecosn), 
where a =} (6+ c), ande = (4 —c)/(6 +); and ~ is equivalent to an integral 
number of circumferences when «= c, and augments by half a circumference 


when tv, next following, attains the value 6. Thus #, like 4, augments continu- 
ously. We have 


6—x=ac(1+cosu), 1+£—c=ae(l —cosy), 
— +) —c)] = ae 


ax = aesinudnu. 


Therefore at = Rdu. 


| 
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As & is a one-valued function of «+ or of a (1 — e cos %), it can be expanded 
in the following periodic series 


R= + a, cos + 2a, cos 2u + 34a,cos3u +... |], 


N, 4, Uy, etc. being constants, the first having the value 


I I 
= f Rau 
oO 


Then ¢ being an arbitrary constant, 
n(¢t+c)=u + a, sinw + a, sin 2” + a, sin . 

This series serves for determining ¢ when + or # is given; but, more fre- 
quently it is + or « which is required in terms of 7 It is necessary, then, to in- 
vert the series. The coefficients of the inverted series are most readily found by 
means of definite integrals. Let us suppose that it is required to find the periodic 
series, in terms of ¢, for a function of + and 1, which we will denote by U. This 
function we assume to be always finite and continuous. The base of hyperbolic 
logarithms being ¢, let us put 

s=e 
and, for brevity, 
25 = a,(s —s~—') + a,(s?—s—*) + 


The equation connecting < and s is 


== 
00 
We can suppose that 
27 27 
Then A, = = Us-*e-# nRdu. 
4 


U being = F(x, 4), we have 
ae sin u | 


U=F — €cos — 


=F 


Supposing that U is reduced to +, it is plain that the coefficient of 2‘, in the 
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development of 1 in powers of 2, is the same as the coefficient of s' in the devel- 
opment of 


in powers of s. 


By adopting the Besselian functions /,"', we have 

i= — j=— 
and the expression, given above, can be written 
(1 
2 } = 
ae|(s+s 


a4, —- — sts 


2 


where 4, = — 

However, unless the coefficients 4, 4, 4... decrease rapidly, this will not 
be a practical method of developing + in a periodic series. Generally it will be 
shorter to employ mechanical quadratures in obtaining the value of the definite 
integral. Let us suppose that 


= $4, + cos + cos + (4,008 37 + . 
Then 


2 
A= 2 fe cos df 
2 
faci —ccosu)[1 + 4, cos + 24, cos + 34, COS +... ] cos du, 


where, to obtain the value of £ corresponding to a given value of ~, we employ 
the equation 
=u + 4,sin + a, sin 2% + 


It will be seen that this method is applicable to a much wider range of ques- 
tions than the motion of planets in elliptic orbits. And the superiority of the 


; 
i! 
| + 
| 
| 
| 
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method of definite integrals over Lagrange’s Theorem for the inversion of the 
series is quite manifest. 
Il. 


In order to illustrate the preceding general theory, let us treat the problem of 
z material points moving about a centre under the action of central forces admit- 
ting a potential which is a function of the sum of the squares of the radi vectores. 
Each point will then move in a fixed plane and its radius vector will describe 
equal areas in equal times. Thus all will be virtually known in reference to 
these motions, provided we are able to express the radi vectorcs as functions of 


the time. 
Let the radi be denoted as 7,, 7,,... 7”,, and the orbit longitudes, measured 


each from any point in its plane, as 4), 4,,...4,. For brevity, put 


Then, if the potential is represented by f(y), we shall have the two equations, 
representing generally all the equations of the problem, 


a*r, ai? 


h,; being the constant of areolar velocity. Consequently if we put 
‘ 


the general form of the differential equations determining the radi vectores will be 


They have the integral, corresponding to that of living forces, 


2 
y 


=2(24 C), 


C being an arbitrary constant. Also we may derive 


dt? er, 


By adding the last two equations, 


dw 


oy of'(n) + 2C 


dh, hh, 
ator” 
at? or, 


an equation involving only the dependent variable ». Multiplying it by the fac- 


ado 
tor 2p 7 , and integrating, we get, 4 being an arbitrary constant, 


at 


To simplify this, we will adopt an auxiliary variable ¢', such that 


150 HILL. ON DIFFERENTIAL EQUATIONS WITH PERIODIC INTEGRALS. 


pea = 2? + 


Inverting this we shall have p as a function of ¢ 
By dividing the penultimate equation by »? and differentiating, we get 


at? p : 


The general equation determining 7; is 


2 2 

a*r, f h; 

at? rs 


As pis now a known function of /, 7; is the only unknown in it, and, conse- 


quently, the equation by itself suffices for determining it. To put the equation 
in a form suitable for integration, let us eliminate f’(7) between the last two 


equations. We get 


d (pdr, — rp) _ 
at? 


A Ado 
— 


a; being the arbitrary constant. By putting 


=) (4%), 


v. 

— 
Whence, by integration, we derive 

fate) 
ay A r2 ’ 
1 
1. 
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au, 


we get a} = 


For convenience, adopting a new constant ¢,, in place of 4;, such that 42/1? 
= a,?(1 — e) the quantity under the radical sign becomes 

[a, (1 u;| [ —a,(I — 
Thus, putting #, = a,(1 -— ¢, cos ¢,), ¢; being a new variable, we get di) = Lids, 
and thus ¢; = 2¢' 4+ 4,, 4; being a constant. Thus we have, in fine, 


=} ja — cos (2¢ + 24;)]}. 


As we have -==I, 


the constants a,, ¢,, and 4, satisfy the relations 
Ja=1, Zae;cos2a,=0, sin 24,0. 
We thus have 2” independent arbitrary constants introduced by integration; the 


number there should be. 
In order to find an expression for the longitudes, we take the general equa- 
tion 
hat 
~ [1 — ¢ cos + a,)] 
— ds 


1—ecos2( + 


ai, 


The integral of which gives 


tan (4, + = \ tan + 


7; being the arbitrary constant. 

To simplify the equations which givé ¢ + ¢ and ¢', we suppose that a(1 + ¢) 
is the maximum value of », and a(1 — cc) its minimum value. Then we can 
adopt a variable ¢ such that 

— cos 
Thus dy = ac sin ¢dz, and we may put 
[f(e) + C] — = sin’ «, 


where RX remains constantly positive throughout the motion of ». Then 


cH fide, 


A 
h — § —— de 
ds. 
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X, being a function of », is also one of a(1 — ¢ cos ¢), and thus is capable of being 
expanded in a converging series of terms, each consisting of a constant multi- 
plied by the cosine of a multiple of ¢. Also » A and 4 / »R can be expanded in 
similar series. Then the period 7, in which » goes through the round of its val- 
ues, is given by the definite integral 

27 


fali—rcose) , 
7 =f ds, 


oO 
and the augmentation of the variable ¢’, in the same time, will be equivalent to 


the definite integral 


J a(t —ecose) 


If the value of the latter is 27, ¢¢ will augment by a circumference while » goes 
through its period. This is the case when f(v) = y; but, in general, this con- 
dition is not fulfilled. 

Provided that <I’ is a positive quantity, it is plain that, after ¢ has gone 
through its period, the longitudes and latitudes, whether as seen from the centre 
or from any of the points, all return to the same values. The same thing is true 
of the ratios of the radi vectores. Thus the movement of the system may be 
conceived as taking place under the operation of two distinct causes. The first 
producing a revolution of all the points about the centre in closed curves and in 
the same time, while the second, having a different period, changes the scale of 
representation of the system in space. 

In the preceding treatment we have supposed that 4? is a positive quantity. 
When this is not the case, some modifications must be made. Let us suppose 
first that =o. Then we have 


do 
» 
+ C]} 


and we may assume 


do 
Sie (2[f(v) + C]} 


| 
\ 
| 
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2 
Also di, = on dy 


tan (A, + = + 


In the second place, let 4? be negative. Here it is only necessary in some 
places to accomplish the integrations by the aid of hyperbolic cosines instead of 
circular. 

The differential equations of this problem, in the case where the radii are 
supposed to describe no areas, were first integrated by Binet.* But the addition, 
to the forces, of the terms arising from centrifugal action, much enhances the 
interest of the problem. 


*See Liouville, Journal de Mathématigues, First Ser. Tome I. p. 457- 


ON THE FOCAL CHORD OF A PARABOLA. 
By Pror. R. H. Graves, Chapel Hill, N. C. 


Let »* = 4ax be the equation to a parabola, s its focus, and PS?’ a focal 
chord. Let the tangent and normal at ?’ meet the diameter through / at .J/and .V. 

It may be easily proved that P17 = P.V = FPP’ and that a similar property 
holds for the tangent and normal at 7’. 

Therefore, if two equal rhombs be constructed on /’?/”’ having two other 
sides of each parallel to the axis, their diagonals are tangents and normals at ? 
and P’; and the tangent at one point is parallel to the normal at the other. 

Each normal chord divides the other in the ratio 1: 3. : 

The chord joining the other ends of the normal chords is parallel to //”’ and 
three times as long. 

A line perpendicular to /P’ at S, and terminated by this parallel chord and 
the pole of PP’, is divided by S in the ratio 1 : 4. 

Hence the locus of the foot of the perpendicular dropped from .S on the par- 
allel chord is a right line, whose equation is + = ga. 

Hence the envelope of the parallel chord is a confocal parabola, having for 
its equation = 32a (ga — +). 

It cuts the original parabola orthogonally where it is cut by its evolute. 
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THE LOGICAL FORM OF GEOMETRICAL THEOREMS. 
By Dr. ALEXANDER MacraRLANe, Austin, Texas. 


In the introduction to the Syllabus of Geometry prepared by the Geomet- 
rical Association of England, there is a portion devoted to the logic of the sub- 
ject; so also in the introduction to Halsted’s Elements of Geometry. The Asso- 
ciation represents the typical theorem by 

If A is B, then Cis D; 
while Halsted takes 
vis 
I have looked at this subject from the point of view of the Algebra of Logic.* 
The typical form which it suggests, has this great advantage over the preceding, 
that we can, by means of it, prove the truth of the Rules of Contraposition and 
Conversion. Let A denote the subject of the theorem (as triangle), x the attri- 
bute of the antecedent, y the attribute of the consequent; then the typical theo- 


rem is ‘ 
2A {x= vy}, (1) 


where » is used to denote a//. 


The contrapositive is 
<A {1 —y =(1 —y) (1 — 2)}. (2) 
Now (2) is true when (1) is true, for it is 
2A {1 
which, when (1) is true, reduces to an identity. 
The converse is 2A{y= 3x}, (3) 
which does not follow from (1), unless 2A {y = x}. 
The obverse is 
2A — 4S (1 —4)(1 (4) 
When the antecedent of the theorem involves two independent attributes, the 
type is 


2A {xy = xyz}. (5) 

This form has two contrapositives ; namely 
YA (1-2) (1 —y)}, (6) 
and YA (1 — — (7) 


both of which are true, when (5) is true. 


* Principles of the A/yebra of Logic, by A. Macfarlane, D. Sc., Edinburgh, 1879. 
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This type has also two converses, and two obverses. 

In the works referred to, a rule of conversion is enunciated, but not proved; 
thus (in Syllabus) “If of the hypotheses of a group of demonstrated theorems it 
can be said that one must be true, and of the conclusions that no two can be 
true at the same time, then the converse of every theorem of the group will nec- 
essarily be true.” 

Let the group of theorems be 


{y= HN), (8) 
2A (10) 
it is also given that X, + % + % = 1; 
and =O, =O, 137, =O. 


By adding (8), (g), (10), we get 

Multiply by ,, then Je 
Similarly MA 
Hence it follows that 


=A = - 
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A SHORT DEMONSTRATION OF THE EXPONENTIAL THEOREM. 
By Mr. J. M. ScHarperte, Ann Arbor, Mich. 


To demonstrate the Exponential Theorem without the aid of the calculus, 
various methods have been devised by mathematicians; but all of the methods 
known to the writer are much longer than the one here given, which, for brevity 
and directness of proof, leaves but little to be desired. 

Assuming that the equation 

(1) 
is true for all values of 1, then ¢, must equal unity. Squaring (1) we have 


a = 1 + 26,4 + (2c, + + + + + 2c, 


(2) 
while, if in equation (1) we substitute 22 in place of x, we have 
= 1 + + + + (3) 
Equating the coefficients of the like powers of x in (2) and (3), we at-once obtain 
c 2 C 3 c 4 


and if the succeeding values of c follow the same law, 


n 


(a — 1)!’ 
which is readily proved by equating the coefficients of the th power of x in the 


second members of (2) and (3). Thus 
+ 26n + 26, + 2c, ~3¢3t int 1) Chin —1) (or Cyn") — 


(4) 


which, if we substitute the assumed values of ¢,, ¢, —\, ¢,—2, ete., becomes 
n n 


When x is even, the last term within the [ ] is (;;,), and when ~ is odd, the last 
term is 2 (,,,"_)). Now, as the value of the expression within the [ ] in equa- 
tion (5) is always equal to 2", it follows that, in general, c, must always equal 


Cc n 
' and hence 


ul? 
2,-2 8,38 
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A TREATISE ON ALGEBRA.* 


This work, as the writers state in their preface, is meant for students who 
have already an elementary knowledge of algebra and geometry, and as a book 
of reference for teachers. Among its many merits may be mentioned the com- 
pleteness of its contents and the exhaustiveness with which each subject is treated. 
Besides the matter usually treated of in works on algebra, it contains some things 
not so common, such as infinitesimals and integration, and brings in at the very 
beginning the idea of function. It also treats of logarithms, continued fractions, 
etc. in a more extended and helpful manner than usual. The use of checks 
wherever available and of methods of shortening work, as in detached coefficients, 
are worthy of note. 

The defects of the work as compared with its merits are small. Among 
them may be mentioned the absence of an index, and the not always clear state- 
ment of a proposition or definition. However, the examples accompanying the 
definitions are very helpful. The work will be of value as a book of reference, 
but as a text-book it is doubtful whether it can be successfully used, unless pre- 
ceded by a complete course of elementary collegiate mathematics. F.P. 1. 


*4A Treatise on Algebra. By Profs. Oliver, Wait, and Jones, of Cornell University. Ithaca, N. Y.: 
Dudley F. Finch. 
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127 | 

A vesseEL of depth a, the top and bottom of which are horizontal planes, is 

filled with a transparent fluid, the refractive index of which at a depth s below 

the surface is 1 + s/a. Two small holes being made in the top, a ray of light 

enters at one hole, is reflected at the bottom, and emerges at the second hole; 
show that the distance between the holes must not be greater than 


2a log, (2 + 4 3). [ Hoover. | 
SOLUTION. 


Let ~ be the declivity of the curved ray at depth «. Then the algebraic ex- 
pression of the hypothesis is 


cos(u+du) 
COs a+2+ ds’ 
which is reducible to tan = 
a+s 


The first integral of this equation is 
log sec w= C + log(a + s), 


and as the range is greatest when y vanishes with 


C= — loga, 
sec 


If in this we put for its value tan~' dz / dx and clear of radicals, we get 


ds 
(a+2)+1[(a+ + 


whence *=alog — 


and the greatest range is found by putting s = @ to be 
24% = 2a log (2 + 7/3). [W. Thornton.| 
[The proposer sends the same solution; except that the greatest range is 
determined by the condition that the ray must not suffer total internal reflec- 
tion on its return to the surface. ] 


128 


A RHUMB line which cuts the meridians at an angle of 30° is projected upon 
a plane tangent to the south pole. The centre of projection is: 1° at the centre 
of the sphere; 2° at the north pole; 3° at an infinite distance on the line of the 
axis. Required the equation of the projection of the curve in each case. 
[ Yale Problems.]} 
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SOLUTION. 


Let 4z — 4 be the latitude of any point on the rhumb line, and ¢ its longi- 
tude, counted from the meridian through the intersection of the curve with the 
equator. Then, if @ be the radius of the sphere, ds an element of the rhumb 
line, and « its constant angle, 

ds = adll sec a 
= asin 4d¢ cosec 4. 


Eliminating ds, dg = tana Sin 0° 
g =tan a log tan }0; 


and the known equation of the rhumb line is 


r=—a 


Let » and ¢ be the polar co-ordinates in the tangent plane of any projected 
point; the pole being at the south pole, and ¢ being equal to the longitude be- 
fore projection, The equations of the projections are, when « = 30°, 


2a 
f 
= 2a tan = 2a V3-%, 
2a 
3° v= asin 


[Frank Miller.] 
144 


Snow that the diurnal path of the shadow of the top of a vertical rod on a 
horizontal plane will be a rectangular hyperbola, if 


2 0 == sin [E+¢] + sin 


where ¢ is the latitude of the place and 0 the declination of the sun, assumed to 
be constant. [R. H. Graves.] 
SOLUTION. 
In the spherical triangle formed by the pole P, the zenith 7, and the position 
of the sun S, we have S? = go° — 0, ZP = 90° — ¢. 
At sunrise and sunset the end of the shadow of the rod will be at infinity, 
i. e. SZ == 90°. 


The asymptotes of a rectangular hyperbola are at right angles, and conse- 
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quently radi vectores to the points at infinity which are parallel to them, also en- 


close a right angle. 
The triangles for sunrise and sunset are similar, and the radz vectores being 
parallel to tangents to S7 at 7, Z must equal 135°. Hence by Napier’s rules 


sin S = cos 45° cos ¢ 
= h[sin(45° + ¢) + sin(45°—¢)]. 


[ A. Lane.| 
EXERCISES. 
151 
INTEGRATE the differential 
sin 4 + cos 4 WW 
sin‘ # + cos*d‘ [E. Frishy.] 


152 
Ir x and / be radii vectores in a parabola, and 2/ the difference of the true 
anomalies, the area of the parabolic sector between the radiz is 


sinf.(r + + cos). [A. Hall.] 


153 
THE normals to a parabola at the points 7, Q, RX pass through the point O 
and meet the axis at the points ?,, Q,, A,; S is the focus and - is the length of 
the /atus rectum; show that 


RR, = SO. [R. H. Graves.) 
154 
Finp the locus of the point from which the four normals drawn to an equi- 
lateral hyperbola form an harmonic pencil. [R. H. Graves.] 
155 


Al TELLS the truth 4 times out of @; 4, whose veracity is unknown, joins A 
in affirming the occurrence of an event. What is the probability that the event 
did occur ? [Artemas Martin. | 

156 

A HAWK which is 100 yards vertically above a hare that is running with a 
uniform velocity in a straight line on a level plane, starts in pursuit, and flying 
continually toward the hare with a uniform velocity, catches the hare at a point 
100 yards distant from where the hare was when the hawk started. What is the 
length of the curve described by the hawk while in pursuit. [/. £. Hendricks.] 
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